HARNACK'S INEQUALITY FOR FRACTIONAL NONLOCAL EQUATIONS 
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\^J • Abstract. Wc prove interior Harnack's inequalities for solutions of fractional nonlocal equations. 

Cn ' Our examples include fractional powers of divergence form elliptic operators with potentials, op- 

(-H , erators arising in classical orthogonal expansions and the radial Laplacian. To get the results we 

r^ ■ use an analytic method based on a generalization of the Caffarelli-Silvestre extension problem, the 

Harnack's inequality for degenerate Schrodinger operators proved by C. E. Gutierrez, and a trans- 
ference method. In this manner we apply local PDE techniques to nonlocal operators. On the way 
a maximum principle and a Liouville theorem for some fractional nonlocal equations arc obtained. 



1. Introduction 



r^ I Very recently, a great deal of attention was given to nonlinear problems involving fractional integro- 

jrt ' differential operators. These problems arise in Physics (fluid dynamics, strange kinetics, anomalous 

transport) and Mathematical Finance (modeling with Levy processes), among many other fields, see 
for instance [5, 6, 8, 15, 21, 22] and the references therein. The main question is the regularity of 
solutions. One of the tools that plays a crucial role in the regularity theory of PDEs is Harnack's 
CS '. inequality, see for example [6, 7, 9, 10, 11, 12, 19, 23, 25, 28]. 

In this paper we show interior Harnack's inequalities for solutions of nonlocal equations given by 
fractional powers of second order partial differential operators. The operators wc consider arc: 
lO , • Divergence form elliptic operators C = — div(a(x)V) + V{x) with bounded measurable coef- 

ficients a{x) and locally bounded nonnegative potentials V{x) defined on bounded domains; 
CO . • Ornstcin-Uhlenbcck operator Ob and harmonic oscillator Hb on M"; 

^D ' • Laguerre operators Lq, L^, L^, LJ^ and L^ on (0, oo)" with a <E (— l,oo)"; 

• Ultrasphcrical operators Lx and Ix on (0, n) with A > 0; 

• Laplacian on domains O C M"; 

• Bcssel operators A a and 5a on (0,C!o) with A > 0. 
For the full description of the operators see Sections 3, 5 and 6. In general, all these operators L are 

?—( i nonnegative, self-adjoint and have a dense domain Dom(L) C L^{fl,dri), where il C R", n > 1, is an 

open set and drj is some positive measure on J7. In Section 2 we show how the fractional powers L'^ , 
< (7 < 1, can be defined by using the spectral theorem. 

Theorem A (Harnack's inequality for fractional equations). Let L be any of the operators listed 
above and < cr < 1. Let O be an open and connected subset of fl and fix a compact subset K <Z O. 
There exists a positive constant C , depending only on a , n, K and the coefficients of L such that 

snpf <C inf /, 

K K 

for all functions f £ Dom(L), f > in fl, such that L'^ f — Q in L^(0,dri). Moreover, f is a 
continuous function in O. 
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Theorem A is new, except for three cases: the Laplacian on K" ([7, Theorem 5.1] and [16, p. 266]), 
the Laplacian on the one-dimensional torus [19, Theorem 6.1] and the harmonic oscillator [23, The- 
orem 1.2]. Harnack's inequality is well-known for divergence form Schrodinger operators with locally 
bounded potentials [12], see also [9, 11, 28]. For the non-divergence form operators listed above the 
result can be obtained by using our transference method of Section 4. Very recently a Harnack's 
inequality for the fractional Laplacian with lower order terms was proved in [25]. 

A novel proof of Harnack's inequality for the fractional Laplacian was given by L. Caffarelli and L. 
Silvcstre by using the extension problem in [7]. Let us briefly explain it here. Consider / : R" -> M. 
as in the hypotheses of Theorem A. Let u{x,y) be the extension of / to the upper half space R"^^ 
obtained by solving 

fdiv(yi"2<Tyy) ^ Q^ ij^ ^n X (0,00); 

\u{x,0) = f{x), onM". 

Let u{x,y) — u{x, \y\), y e R, be the reflection of u to R"+^. The hypothesis (— A)'^/ = in O implies 
that y^~'^'^Uy{x,y) — >■ as y — J- O"*", for all x G O. This is used to show that u is a weak solution of 
the degenerate elliptic equation with A2 weight 

div(|j/|l-2<Ty ~) ^ Q^ ij^ X (_^^ ^) ^ jjn+1^ 

for some R > 0. Recall that a nonnegative function u on M" is an A2 weight if 

1 /• \ / 1 , _i 



bZ\w\Jb1[w\' ^ '<°^- 

Then the theory of degenerate elliptic equations by E. Fabes, C. Kenig and R. Scrapioni in [10] says 
that u satisfies an interior Harnack's inequality and it is locally Holder continuous, thus f{x) ~ u{x, 0) 
has the same properties. 

The idea of [7] was also exploited in [23] for the case of the fractional harmonic oscillator (— A-f |a;|^)'^ 
on R", under the additional assumption f G C^. In [23] a generalization of the extension problem 
was proved that applies to a general class of differential operators, and it was used to get the result 
for the harmonic oscillator. We observe that, instead of the theory of [10], Harnack's inequality for 
degenerate Schrodinger operators of C. E. Gutierrez [12] had to be applied. 

To get Harnack's inequalities for fractional powers of the operators listed above we push further 
the Caffarelli-Silvestre ideas. We proceed in two steps. First we use two tools: the extension problem 
of [23] and Harnack's inequality for degenerate Schrodinger operators of C. E. Gutierrez [12]. These 
are enough to get Theorem 3.2, from which the result for divergence form elliptic operators with 
potentials and some Schrodinger operators from orthogonal expansions is deduced. Secondly, we 
apply systematically a transference method that permits us to derive the results for other operators 
involving terms of order one and in non-divergence form. The transference method is inspired in ideas 
from Harmonic Analysis of orthogonal expansions, where it is used to transfer L^ boundedness of 
operators, see for example [1, 2, 14]. In that case, the dimension, the underlying measure and the 
parameters that define the operators play a significant role. Here we can obtain our estimates without 
any restrictions on dimensions or parameters. 

Let us remark that in Theorem A we require the condition / > all over Q, which is needed to 
ensure that the solution to the extension problem u is nonnegative in i7 x (0,oo). In fact, u can 
be given in terms of the solution e~*^/ of the L-heat diffusion equation, see Theorem 2.2 below, so 
we only would need the condition e~*^/ > in C Certainly it is sufficient to assume that e~*^ is 
positivity-preserving (see (2.2) below), but this hypothesis is not strictly necessary. 

As a by-product of our method, we obtain a Liouville theorem for fractional powers of divergence 
form elliptic operators on R", see Remark 3.3. We also get a maximum and comparison principle for 
general fractional operators, see Remark 2.1. 

In Section 2 we present the definition of fractional powers of differential operators, we get maximum 
and comparison principles and we state the extension problem of [23]. The method of reflections for 
proving Harnack's inequality for divergence form elliptic Schrodinger operators is given in Section 3. 
The transference method is explained in Section 4. 
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The rest of the paper is concerned with the proof of Theorem A in each case. As the reader 
may notice, we have two sets of apphcations of our method: operators with discrete spectrum and 
operators with continuous spectrum. In the first set we have divergence form elhptic operators in 
bounded domains and classical operators related to orthogonal expansions in possibly unbounded 
domains (Sections 3 and 5). In the second set (Section 6) we have the Laplacian (Fourier transform) 
and the Bessel operator (Hankel transform), that generalizes the radial Laplacian. 

We will present most of the results about Harnack's inequalities in the case when the sets K and 
O in Theorem A are balls inside fi. In that situation the constant C does not depend on the radius 
of the balls. By the standard covering argument [11, Theorem 2.5] the general result can be easily 
deduced. 

Through this paper we always take < cr < 1. 

2. Fractional operators and extension problem 
Along this paper all the operators will verify the following 

General assumption. By L = L^ we denote a nonnegative self-adjoint second order partial differ- 
ential operator with dense domain Dom(i) C L^{Q,dr]) = L^(J7). Here il is an open subset of M." , 
n > 1, and drj is a positive measure on Q. The operator L acts in the variables x G M". 

The Spectral Theorem can be applied to an operator L as in the general assumption, see [20, 
Chapter 13]. Given a real measurable function (j) on [0,oo), the operator (j){L) is defined as (/)(L) = 
L (/)(A) dE{X), where E is the unique resolution of the identity of L. The domain Dom(0(L)) of 0(i) 
is the set of functions / e L'^{n) such that J„ |^(A)| dEff{X) < oo. 

In this paper we are going to use: 

• The heat-diffusion semigroup generated by L, defined as <^(i) = e~*^, i > 0. For / S L^{il), 
we have that v ~ e^*^/ solves the evolution equation Vt ~ ~Lv, for t > 0. Moreover, 
l|e-*^/||L2(n) < ||/||L2(n), for all t > 0, and e"*^/ ^ / in L^{n) as t ^ 0+. 

• The fractional powers of L, given by 4>{L) = L°' , with domain Dom(ij'^) D Dom(L). When 
/ e BomiL") we have L^e'^^f = e-*^L"f. If / G Dom(L) then {Lf,f) = [[Li/^j 
where (•, •) denotes the inner product in L^(J7). Also, for / e Dom(i), 



iL2(n), 



1 C°° dt 

(2.1) L^f{x)^^,^J^ (e-*^/(^)_/(a:))^^, h, L^n), 

where F is the Gamma function, see for example [29, p. 260]. 
We will usually assume that the heat-diffusion semigroup e~*^ is positivity-preserving, that is, 

(2.2) f >Oonn implies e"*^/ > on O, for all t > 0. 

Remark 2.1 (Maximum and comparison principle for £'^). Let L be as in the general assumption. 
Under the additional hypothesis (2.2), the following comparison principle holds, li f,g G Dom(i), 
f > g in D, and /(xq) = g{xo) at a point xq G fl, then L"^ f{xo) < L'^ g{xo). This comparison principle 
is a direct consequence of the maximum principle: if / G Dom(L), / > 0, /(a;o) — 0, then L'^f{xo) < 
(for the proof just observe in (2.1) that T{—a) < and e^*^/(a:o) > 0). 

Theorem 2.2 (Extension problem [23, Theorem 1.1]). Let L be as in the general assumptions and 
f e Dom(L'^). Let u be defined as 

, X V^" r tLr, . "' dt 



...^ 4-F(a)7o ''"' t^+" 

' e-*^(LV)(^)e-« 



T{a)Jo ''''-' t^ 
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for a; e f2, 2/ > 0. Then u G C°°((0, oo) : Doni(L)) n C([0, oo) : L^(J7)) and it satisfies the extension 
problem 

J-L^U+^^Uy+Uyy=0, X G fi , 7/ > , 



u{x,0) = f{x), xefl. 



T-1/2 



In addition, for Ca = r(i-a-) ^ ^' 

(2.5) -c, lim yi-2a^.y(x,y) = LV(a:^). 

y->0+ 

We must clarify in which sense the identities in Theorem 2.2 are taken. The first equahty in (2.3) 
means that for any g £ L^(ri), 

(-(•,2/),5(-)) = ^/ (e-*^/,5)e-«r_, , > Q, 

and similarly for the second one. Also (2.4) in general means that {^—^Uy{-,y) + Uyy(-,y),g{-)) = 
{Lu{-,y),g{-)), for aU y > 0, with {u{-,y),g{-)) — >■ (/,<?), as y — )■ 0+, and analogously for (2.5). By the 
second identity of (2.3), a change of variables and dominated convergence, we have 

limsup||yi-2-^,(z,y)||i.(^,)<— — limsup / We-^^L^ f)\\hin'f''' Z^ 



,-ii 



c„ 



L" f\\L^{n')7 for a-ny measurable set i7' C 51. 



3. Harnack's inequality for fractional Schrodinger operators 

In this section we consider a uniformly elliptic Schrodinger operator of the form 

£ = - div(a(a;)V) + V, on fl C W\ 

Here a = (a*^) is a symmetric matrix of real-valued measurable coefficients such that /i^^'^l^fp < 
a,{x)^ ■ ^ < fJ.\^\'^, for some constant fi > 0, for almost every x G 51 and for all ^ e W\ The potential 
y is a locally bounded function on 51. Here 51 can be an unbounded set. We assume that C satisfies 
the general assumption at the beginning of Section 2, with drj{x) = dx, the Lebesgue measure. The 
domain of C is Dom(£) = Wq' (51) n i^(52, V{x) dx). The Sobolev space Wq' (fi) is the completion 
oiC^ifl) under the norm ||/||^i,2(o) = il/lli2(n) + W^ f\\h(ny ^0*° t^at Dom(/:) is dense in L'^{n). 
For / e Dom(£), 



{Cf,g)^ / {a{x)Vf-Vg + V{x)fg)dx, g e W^'\n) n L'{n,V{x) dx). 
Jn 

Theorem 3.1 (Reflection extension). Fix a ball Bpi{xa) C 51, .tq G 51, _R > 0. Let it : 51 x [0, i?) — > K 
be a solution of the extension equation in (2.4) with L = C in Bii{xo) x (0, i?). Define the reflection 
of u to ft X (—R,R) by u{x,y) = u{x, \y\), x Gft, y G {—R,R). Suppose that 

(I) lim,y^o+ \\y^^'^'^Uy{x,y)\\L2(^BR{xo),dx) =0; and 
(H) \\'^xu{x,y)\\L2(Bii{xo},dx) remains bounded as y -^ 0+. 
Then u verifies the degenerate Schrodinger equation 

(3.1) div(|y|i-2'^6(a;)Vw) - \y\^~^''V{x)u = 0, 

in the weak sense in B := Ux,y) G R"+^ : |x — .Top + y^ < R^} , where the matrix of coefficients 
b = {b^i) is given by V^ = a'^ 6"+i-J = &*'"+i = 0, 1 < i, j < n, and 5"+i^"+i = 1. 

Proof. Let (f G C^{B). Take any < S < R. Since w is a solution of the extension equation in (2.4) 
for £, for any fixed y G {S, R), we have 

/ {aix)VxU-Vx^ + Vix)uip)dx= [ \yf''-^dy{\y\^-^''uy)^dx. 

JBr(xo) JBaixo) 
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Recall that we are assuming that u £ C°°((0,i?) : Dom(£)). By integrating the last identity in y, 
applying Fubini's theorem and integration by parts, 

R . 

Ij/I^"^"^ / [a[x)\/ xU ■ \/ x'f + V [x)uLp) dx dy 

s Jbr{xq) 

6^^^''uy{x,6)ip{x,6)dx- / \y\^~'^'^Uy{x,y)ipy{x,y)dydx. 

r{xo) J Bftixo) JS 

From here we get 



(3.2) / {bix)\7u-\/ip + Vix)uip)\y\^'^''dxdy 

JBR(xo)x{\y\>S} 



6^ ^''uy{x,S)ip{x,~S)dx ~ / S^ ■^''uy{x,S)ip{x,S)dx. 

r{xo) JBr{xo) 

We are ready to prove that m is a weak solution of (3.1) in B. We have to check that 



/:= / {b{x)Vil-Vip + V{x)ilip)\y\'-^'' dxdy = 0. 
Jb 



By using (3.2), 



I=\ I + dxdy 

3n{\y\>S} JBn{\y\<5} J 

S ^ '^Uy[x,5)ip{x,~5)dx ~ I 5^ '^Uy{x,S)Lp{x,S) dx 



Br{xo) •'Br{xo) 

V{x)uip\y\^-^'^dxdy. 



+ [ b'^Vn-Vip\y\^-'^'^dxdy+ I 

JBn{\y\<S} Jsn{|y|«5} 

As (5 — > 0+, the first and second terms above tend to zero because of (I). Also the fourth term goes 
to zero because V [x)u\y\^^'^'^ € L\^^. Since \\^ xu{x,y)\\j^2ig/^ \ ^^\ remains bounded as y -> 0+, for 
any small (5 > there exists a constant c > such that if \y\ < S then ||Va:u(a;, y)||^2f3 ,^ i ^^l < c. 
This property and (I) imply that the third term above tends to zero as i5 — >■ 0"''. D 

Theorem 3.2 (Harnack's inequality for C^). Let C be as above. Assume that the heat- diffusion 
semigroup e~*^ is positivity-preserving, see (2.2). Let f G Dom(£) be a nonnegative function such 
that L" f ^ in L"^ [B nixo) , dx) for some ball Bii{xo) C fl. Suppose that \\'Vxu{x,y)\\L2(BR{xo}Mx) 
remains bounded as y — > 0'*', where u is a solution to the extension problem (2.4) for C and f. There 
exist constants Rq < R and C depending only on n, a, /i, and V, but not on f, such that, 

sup f <C inf /, 

Br ^^ 

for any ball Br with B^r C Bii(xq) and < r < Rq. Moreover, f is continuous in Bji{xif). 
In order to prove Theorem 3.2 we use Theorem 3.1 and the following version of 

Gutierrez's Harnack inequality for degenerate Schrodinger equations. Consider a degener- 
ate Schrodinger equation of the form 



(3.3) -6iY{d{X)Vv) + V{X)v = id, Xe 



bN 



where a = (a'-') is an A^ x A^ symmetric matrix of real- valued measurable coefficients such that 
A-^a;(X)|^P < diX)^ ■ ^ < \ujiX)\^\^, for some A > 0, for almost every AT € M^ and for all ^ S M^. 
The function lj is an A2 weight. The potential V satisfies V/cu G L£ locally, for some large p ~ pn,ui- 
Let O be any open bounded subset of M.^ . Then there exist positive constants r'o,7 depending only 
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on A, A^, w, O and V such that if v is any nonnegative weak solution of (3.3) in O then for every ball 
Br with Bsr C O and < r < ro we have 

sup V < ^ inf V. 

Br/2 ^-/^ 

As a consequence, v is continuous in O. See [12]. 

Proof of Theorem 3.2. Since C^ f = in L^{Bii{xo),dx), by (2.6) and the hypothesis on VxU, we 
see that u satisfies the conditions of Theorem 3.1. Now, equation (3.1) is a degenerate Schrodinger 
equation with A2 weight uj{x, y) — |y|^~^'^ and potential V = \y\^~'^'^V{x) such that V /lo £ L^ locally 
for all p sufficiently large. By Gutierrez's result just explained above, Harnack's inequality for u holds. 
By restricting u to y ~ Q we get Harnack's inequality for /. Moreover, u is continuous in Bji{xo) and 
thus /. D 

3.1. The case of nonnegative potentials. Under the additional assumptions that fi is a bounded 
set and that the potential V^ is a nonnegative function in J7, we can prove Theorem A for C^ . In this 
case the domain of C is Dom(£) = W^ ' (17) and it is known that e^*^ is positivity-preserving, see [9, 
Chapter 1]. Let / G Wo'^(17), / > 0, such that £'"/ = in L^{BR{xo),dx) for some ball Br{xo) C Vt, 
R > 0. Denote by u the solution of the extension problem for / as in Theorem 2.2. By virtue of 
Theorem 3.2, to prove Harnack's inequality for £°' we just have to verify that u satisfies condition 
(II) of Theorem 3.1. As / e W^'^^n), by the ellipticity condition, 

(3.4) A^-'i|V/||i.(^,,,) < / a{x)Vf-Vfdx<{CfJ) = \\C'^' .f\\l^n,d.), 

(for the last equality see Section 2). Now, since u g C^((0, 00) : Wq' {^)), S/xu{x,y) is well defined 
and belongs to L'^{VL,dx) for each y > Q. We can apply (2.3), (3.4) and the properties of the heat- 
diffusion semigroup e~*^ stated at the beginning of Section 2 to get 

,,2<T /-oo 2 7, 



\\'^xu{x,y)\\L2i^Bn(x«).dx)<-r-fi—: l|Ve ''f\\L^n,dx)e « 



< 



< 



4-r(CT)yo " ^Ui.~,..,axr ^1+a 

dt 



1/2 y / 11 -i£/'l/2 /-ii 



4-r(f7)7o JUi.-(^^Mxr ^1+a 



r(a) 



f (1) e-^f = //^||/:V2;||^,^^^^^^. 



Thus \\\'xu{x,y)\\]^2i^Sj^^ri;o),dx) remains bounded as y — >■ 0+ and (II) in Theorem 3.1 is valid. Hence 
Theorem A is proved for this case. Observe that, in particular. Theorem A is valid for the Laplacian 
in bounded domains with Dirichlet boundary conditions. 

Remark 3.3 (Liouville theorem for fractional divergence form elliptic operators). Let fl = M" and 
V = 0, that is, £ = - div(a(a:)V). Take / g Dom(£) = W'^^'^{W). The following Liouville theorem is 
true: If / > on R" and £°'/ = in L^(]R") then / must be a constant function. Indeed, for this /, 
the reflection m of u is a nonnegative weak solution of (3.1) with V^ = in M"+^, so u is constant and 
therefore / is a constant function. Here we have applied the Liouville theorem for degenerate elliptic 
equations in divergence form with A2 weights, which is a simple consequence of Harnack's inequality 
of [10]. 

Remark 3.4. Since our method is based on Gutierrez's result [12], we are not able to get the exact 
dependence on a of the constant C in Harnack's inequality of Theorem 3.2. 

4. Transference method for Harnack's inequality 

In this section we assume that L satisfies the general assumptions of Section 2. We explain in 
detail a general method to transfer Harnack's inequality from L"' to another operator L"' related to 
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L. This method will be useful when considering differential operators arising in classical orthogonal 
expansions and also for the Bessel operator. 

Firstly, by a change of measure, we have the following trivial result. 

Lemma 4.1. Let M(x) € C'°°(fl) be a positive function. Define the isometry operator U from 
L'^{fl,M{x)'^dr]{x)) into L^{fl,dri{x)) as {U f){x) = M{x)f{x). Then if {^Pkjken"^ ** o.n orthonormal 
system in L^{^, M{x)'^dri{x)) then {U(pk}ke'M" *-' also an orthonormal system in L^{fl,dri{x)). 

Next we set up the notation for the change of variables. 

Definition 4.2 (Change of variables). Let /i : fJ — > fj C R" be a one-to-one C°° transformation on 
n. Denote the Jacobian of the inverse map h~^ : f2 — )• O by | J/i-i |. We define the change of variables 
operator W from L'^{fl,M{h~^{x))'^ \Jfi-i\dr]{x)) into L'^{fl,M{x)'^dr]{x)), where dr]{x) — ri{x) dx for 
some positive density 77, as 

iwf)ix)^f{h{x)), x&n. 

Now we are in position to describe the transference method. By using the definition above and 
Lemma 4.1 we construct a new differential operator. This new operator will be nonnegative and 
self-adjoint in L?{p.,dfi{x)), where (l = h{Q) and dfj{x) := M{h~-^{x))^ \Jh-^\ dri{x). Let 

L:= {U oW)-^ oLo {U oW). 

If E is the resolution of the identity of L then the resolution of the identity E oi {U o W) o L verifies 

dEf,g{X) == dE(^uoW)f,{UoW)gW, f,9 & L^{Q,,df]). 

Therefore if / G Dom(Z'^) then we see that the fractional powers of L satisfy 

L^f = {Uo W)-^ oL" o{Uo W)f. 

Lemma 4.3 (Transference method). // Theorem A for U' is true, then the analogous statement for 
L'^ is also true. 

Proof. Let / G Dom(Z'^), / > 0, such that Z'^/ = in L'^{0, dfj), for some open set O C Cl. Take a 
compact set K C O. We want to see that there is a constant C depending on K and L'^ such that 

(4.1) sup/<Cinf/. 

K i< 

Observe that, by the definition of dfj and since dri{x) ~ rjix) dx, 

\L"o{UoW)fix)\^diix)= f \L^f{x)\^dfj{x)=0, 
/i-i(O) Jo 

and ([/ o W)f e Dom(L) is nonnegative. By the assumption on L'^ , there exists C depending on 
h^^{K) and L'^ such that 

sup {U o W)f < C inf ([/ o W)f, 

and {U o W)f is continuous. In particular, / is continuous. Since M{x) is positive, continuous and 
bounded in h^^{K), 

sup Wf < C inf Wf. 

This in turn implies (4.1) as desired. D 

5. Classical orthogonal expansions 

In this section we consider operators L (as in the general assumptions of Section 2) for which there 
exists a family {ipk^ks^^-n. of eigcnfunctions of L, with associated nonnegative eigenvalues {Afe}ji.gj^n, 
namely, Lipkix) = Xkfkix), such that {(pk} is an orthonormal basis of L^ijl, drf). In all our examples, 
the eigenvalues will satisfy the following: there exists a constant c > 1 such that Afc ~ |fc|'^, for any 
k = (fci, . . . , fc„) € Ng , |fc| = fci + • • • + fc„. We also suppose that the eigcnfunctions ipt are in C^(J7) 
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and that their derivatives satisfy the following local estimate. For any compact subset iiT C 51 and 
any multi-index (3 E Wq, \f3\ < 2, there exist e = ex.js > and a constant C ~ Ck,p such that 

(5.1) \\DP^k\\L^(K,d,,)<C\k\\ 



for any k e Wq. For / e L'^{Vt,drj) the heat-diffusion semigroup can be written as e *^f{x) 

^|fel=0 



J2m=o^ *^^Ck^k{x). For < (7 < 1, the domain of L'^ is given as Dom(L'^) = {/ e L^{Q,dri) : 



X]|fe|=o '^ft'^l'-fel^ < °°}' where Ck denotes the Fourier coefficient of / in the basis (fk- Ck = {f^fk) = 
Jq ffk d-q- Given / G Dom(L'^) we have L'^f{x) = J2^\=o K^^k'fkix). 

Under these assumptions we can show that the solution u of the extension problem is classical. To 
this end, let K be any compact subset of fl. First we show that the series that defines e~*^f{x) is 
uniformly convergent in Kx (0, T), for every T > 0. Indeed, by applying that A^ ~ \k\'^, estimate (5.1), 
the inequality s''e~^ < Cpe~^^'^ (valid for s,p > and some constant Cp > 0) and Cauchy-Schwartz's 
inequality, 

= (x)|<-g^^(t^/^|fcne-'^*l'=l' 



L^{n,dij} 




L^{n,d7i) ' 



X e K, 



and the uniform convergence follows. As a consequence, u in (2.3) is well defined, for by the estimate 
above, for any x G K and y > 0, 



°° _tL,, . _Z, dt f^ e-4r dt 



\^-'-f(-)^^^\jTT^ < C ll/IL^(o,..) / y TT^mTT ^^ < Fiy), 



t- 



for some function F = F{y). This estimate also implies that in the first identity of (2.3) we can 
interchange the integration in t with the summation that defines e~*^/(x) to get 



^ '^^'^> \k\>Q -^0 ^ 

By using (5.1) and the same arguments as above, it is easy to see that this series defines a function 
in C^(J7) n C^(0,(X)). Moreover, since each term of the series in (5.2) satisfies equation (2.4) in the 
classical sense, we readily see that u is a classical solution to (2.4). 

Next we will present the concrete applications. 

We will take advantage of well-known formulas, see for instance [1, 2], to apply our transference 
method to get Harnack's inequality for operators of classical orthogonal expansions which are not of 
the form considered in Section 3. A remarkable advantage of the transference method is that we do 
not need to check that the semigroup e~*^ is positivity-preserving. 

5.1. Ornstein-Uhlenbeck operator and harmonic oscillator. In [13], Gutierrez dealt with the 
Ornstein-Uhlenbeck operator 

Ob = -A + 2B2;- V, 

where B is an n x n positive definite symmetric matrix. The operator Ob is positive and symmetric 
in L^(R",d7B(a;)), where d'y^^x) = (det B)"/^7r^"/^e~^^'^da; is the B-Gaussian measure. Let us 
consider the eigenvalue problem Obw = Ait), with boundary conditions w{x) = 0{\x\''), for some 
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fc > as |a;| — > oo. Firstly, let us assume that the matrix B is diagonal, which means that 



B = D 



(di ••• 0\ 
da ■•• 



\{) ••• dn) 



with di > Q for \ < i < n. It is not difficult to see that in this case the eigenfunctions w are the 
multidimensional Hermite polynomials defined by Hj?{x) = Hk^{\/dlxi) ■ ■ ■ Hk^^^fd^Xn), A: G Ng, 
with eigenvalues 2{k ■ d), d = (di, . . . ,dn), where i/fc. is the one-dimensional Hermite polynomial of 
degree ki, see [13]. For the general case, since B is a positive definite symmetric matrix, there exists an 
orthogonal matrix A such that ABA* = D, where A* is the transpose of A. Then the eigenfunctions 
become H^{x) = H^{Ax). 

Let us also consider the harmonic oscillator 

'Hd = -^+\Dx\'^, 

where D is a matrix as above, with zero boundary condition at infinity. Under these assumptions 
Hd is positive and symmetric in L^(R",da;). It is well known that the multidimensional Her- 
mite functions h^{x) = (det £')"/'''7r~"/^e ^Hj^{x), are the eigenfunctions of T-Ld and T-Loh^ = 
{2{k ■ d) + X)"=i ^i) ^k ■ The Hermite functions form an orthonormal basis of L^(R", dx). 
Observe that we may also consider 



hd - y^d,, 

2 = 1 

since it has the same eigenfunctions as Hd with eigenvalues 2{k ■ d) > 0. We can also put a more 
general matrix B in the place of D] we will prove Harnack's inequality for it by using the transference 
method. 

5.1.1. Proof of Harnack's inequality for (Hd)" ■ To show Harnack's inequality for {HdY ^^ have to 
check that all the conditions of Theorem 3.2 hold. 

The potential here is V{x) = |-Da;p, which is a locally bounded function on R". 

By Mehler's formula [13, 24, 26], e~*^° is positivity-preserving. 

In [26], it is shown that there exists C such that ||/i^||L°°(R",dx) ^ C for all k. Using the relation 

2d,^h^{x) = ^,[{2hy/'h^_^^{x) - (2fc, + 2)V2/,f^^^(^)) , 

where e^ is the i-th coordinate vector in Ng , we see that (5.1) is valid for hj^{x). Therefore the solution 
u to the extension problem given in (2.3) for Hd is a classical solution. 

Let / G Dom(HD), / > 0, such that {HdY f = in L'^{BR{xo),dx). We have to verify that 
\\Vxu{x,y)\\^2(^Bii{xo)Ax) remains bounded aay^ 0+. In fact, we wiU have \\yxuix,0)\\^2(^Bn{xo).dx) = 
ll^2:/(2;)llL2(s^(i;o),da;)- Indeed, as we can write / = J2'\k\=o '^khj^ , by (5.2) and the identity for the 
derivatives of the Hermite functions hj^ given above, 

(dxi + VdiXi) {u{x, y) - f{x)) 
^'■'^ -E^^-V^(2fcO^^^/^f-.(.)(4^/°°e-(-+^r...Oe-^^-l). 

Observe that the term in parenthesis above is uniformly bounded in y and, since 

2a /-oo 2 rl+ 1 r°° 2 rim 



4'^r(cr) Jo t^^" r(CT) Jo w 
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we readily see that it converges to when y — > 0+. Moreover, as / e Dom('H_D) 

1/2 

2 y^ c^kidi 1 < oo. 



J2 ckVd,{2hy/^h^_^jx) 

|fc|=0 



^dx) \ 1^1=0 



Hence, by dominated convergence in (5.3), we get that [dxi + \/diXij u{x,y) — >■ (^dxi + ^/dixij f{x) 
in L^(IR",(ix) as y ^- 0+. Since u{x,y) — >■ /(x) in L^(R",(ia;) and \fdiXi is a bounded function in 
Br{xo), we have that \/dixiu{x,y) converges to \/dixif{x) in L^ (B h^xq) , dx) as y — >■ 0"^. Hence 
\7xu{x, y) -^ Va;/(x), as 2/ -> 0+, in L^{Br{xo), dx). 

5.1.2. Proof of Harnack's inequality for (O I)] . We apply the transference method explained in Sec- 
tion 4. For this case we take M{x) = (detD)"/'*7r^"/^e r^ and h{x) ~ x. Clearly h^{x) = 
{U o W)Hj^{x) and we have the relation 

(5.4) OdH^ = {Uo W)-^ o j -Hd - ^ d, J o (t/ o W)H^. 

Sec also [2]. It can be easily checked, as done for {T-LdY above, that the operator {T-Ld — X]r=i ^0 
satisfies Harnack's inequality. Hence the conclusion for {OdY follows from Lemma 4.3. 

5.1.3. Proof of Harnack's inequality for [O-b) ■ Consider the change of variables /i(a;) = A*x and call 
W the corresponding operator as in Definition 4.2. Then it is easy to check that 

OB(Hf o h-'){h{x)) = OoiJf (x). 

Then we have Ob = VF~^ o Ojj o W and the result follows by the transference method. 

5.1.4. Proof of Harnack's inequality for (He) • We observe that parallel to the case of the operator 
Ob we can get T-Lb = W^^ o T-Ld ° W with W as in Subsection 5.1.3 above and then we get Harnack's 
inequality for the operator {T-L-bY ■ 

5.2. Laguerre operators. We suggest the reader to check [1, 14, 17, 24, 26] for the proof of the basics 
about Laguerre expansions we use here. Let us consider the system of multidimensional Laguerre 
polynomials L'^{x), where k £ Nq, a = (ai,--- ,a„) E (— l,oo)" and x £ (0, oo)". It is well known 
that the Laguerre polynomials form a complete orthogonal system in L^((0,oo)", (i7Q,(x)), where 
d'jaix) = x"^e~^i dxi ■ ■ ■ x""e~^" dxn- We denote by Zj? the orthonormalized Laguerre polynomials. 
The polynomials L^ are eigenfunctions of the Laguerre differential operator 

namely, Lq.(LJ?) = 1^1-^/^. There are several systems of Laguerre functions. We first prove Harnack's 
inequality for the operator L^ (related to the system <^^ below) and then we apply the transference 
method of Section 4 to get the result for the remaining systems. 

5.2.1. Laguerre functions (p^. This multidimensional system in L^((0,c»)",d/io(x)), where d^o{x) = 
dxi ■ ■ ■ dxn, is given as a tensor product ^"^{x) = (p^^(x\) ■ ■ ■ (/3^"(xn), where each factor (/)^'(xi) = 
x"'(2xi)-'^/^e~^>/^L^'(x|). The functions Lp% are eigenfunctions of the differential operator 

1 " 1 / 1 

(5.5) LS = J(-A+N^)+Ei^^"' 

Z — 1 * 

namely, 

(5.6) \.l^l{.x) = E (fc. + ^) ^{x^. 
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Clearly, the functions Lp'^ are locally bounded in (0,oo)". Observe that 

(5.7) d..vi{x) = -ifcr/Vfcitr(^) - (^' - ^ ("^ + 0) ^^(^)- 

Therefore, (5.1) holds for this system and we get that the solution u in (2.3) of the extension problem 
for L^ is classical. Moreover, it can be easily seen from [24, p. 102] that e"*'-'" is positivity-preserving. 

Let us prove Theorem A for (L^)"^. We can do this as we did for {'Hd)'^ above by following the 
reasoning line by line, but with some modifications as follows. Let / € Dom(L^), / > 0, such that 
(L^)"^/ = in L^{Bii{xo),dfj,o{x)), and let u be the corresponding solution to the extension prob- 
lem. By (5.7) and a similar argument for that of H"^ we can check that \\'^xu{x,y)\\j^2/g^/^ \ ^^ i^\\ 
converges to \\^xf\\L2iB^(x ) duoix))' as y ^- 0+. Moreover, the potential in (5.5) is locally bounded. 
Hence, by Theorem 3.2, / satisfies Harnack's inequality and it is continuous. 

Note that the same arguments above can be used for (L^ — ^^)°' instead of (L^)*^, so it also 
satisfies Theorem A. 

5.2.2. Laguerre functions i'^. The Laguerre functions i^ are defined as i^ix) = i'^{x\) ■ ■ ■ £^"(a;„), 
where £^' are the one-dimensional Laguerre functions P^^\xi) = e~^'/^L^'(xi). Each i^, is an eigen- 
function of the differential operator 

^'-i:(-.^- (".+<+? 

More explicitly, L^i"^ = X)r=i {^i + ^^^2"'") ^V- ■ ^or d^a{x) = x"^ ■ ■ -x^^dx, the operator L^ is pos- 
itive and symmetric in L^((0, oo)", d/iQ(a;)). The system {£^ : k G Nq} is an orthonormal basis of 
L2((0,oo)",d/^a(x)). 

To apply the transference method we set M{x) = 2"/^x"^ • • • a:"" and h{x) ~ {xf, . . . , .t^J. 

Then t/o VK is an isometry from L^((0, oo)",d^Q(a;)) into L^((0, oo)", d^o(a^)) and L^ = {U oW)^^ o 
L'^o{UoW),scc [1]. 

5.2.3. Laguerre functions tp'^. Consider the Laguerre system 'ip'^{x) = 'ip'^^{xi) ■ ■ ■ip'^"(xn), which is 
orthonormal in L^((0, oo)", d/i2Q+i(a^)), where dfi2a+iix) = x\°'^^'^dxi ■ ■ ■ x'^"'^^dxn and V'^' is the 
one-dimensional Laguerre function V'fe'(2^i) = n/2 ^^'(xf ). The functions "0^ are eigenfunctions of the 
operator 

4 ^ ' ' ^ ^^^ 4xi dxi 

2—1 

In fact, l^ti^t) = Zto {h + ^) C- 

For the transference method we have to take M(x) = x"^ • • • a;"" and h{x) ~ x. Then UoW 

is anisometry from L'^{{0, cx))", d^2a+i(a;)) into L^((Q, oo)", d/io(a;)) and L^ = {UoWy^oL'^o{UoW), 
see [1]. 

5.2.4. Laguerre functions C^. The functions C'^{x) = C'^^{xi) ■ ■ ■ C'^"{xn) form an orthonormal sys- 
tem in L^((0, oo)",d/io(a:)), where £^' is the one-dimensional Laguerre function given by C^Xxi) = 
x" ^^'(si). The functions C^ are eigenfunctions of the operator 



E 



d d Xi 



^ dxj dxi 4 4x, 

i— 1 



In fact, l.i[Ct) = Er=o {k. + ^) Cf^. 

Apply the transference method with M{x) ~ 2"/^a;j^' • • • xj and h{x) ~ (xf, . . . , x-^J. Then UoW 
is an isometry from L^((0,oo)",d/io(x)) into itself and Lf = ([/ o W)^^ o L^ o ([/ o W), see [1]. 
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5.2.5. Laguerre polynomials L'^ . Finally consider the Laguerre polynomials operator L^ . Let M(x) = 
2"/2e-kl /22;"i+ / ■ ■ ■ Xn"'^ and h{x) — (a;^, . . . ,x^J. We have that the operator U o W is a.n 
isometry from L2((o,oo)",d7a(a;)) into L2((0,oo)",d/io(x)) and L„ = {UoW)-^o(L^^^)o{UoW), 
see [1], so the transference method applies. 

5.3. Ultraspherical operators. Here we restrict ourselves to one-dimensional expansions. We de- 
note the ultraspherical polynomials of type A > and degree k G No by P^ix), x e (—1,1), see 
[17, 18, 24]. It is well-known that the set of trigonometric polynomials {P^{cos9) : 6 £ (0,7r)} forms 
an orthogonal basis of i^((0, n),dmx{9)), where dmx{6) — sin^ 9 d9. The polynomials P^(cos 9) are 
eigenfunctions of the ultraspherical operator 

L>^--^.-^>^-ot9±^X\ 

that is, L\P^{cos9) = (fc + A)^P^(cos6'). We denote by P^{cos9) the orthonormalized polynomials 
given by 21--2I Jw^ i 2a) ^ki^^^ ^)- There exists a constant A such that |P^(cos0)| < Ak'^^^^, see 
[18]. This and Stirling's formula for the Gamma function [17] imply that there exists C such that 
|Pfc (^)l < Ck for ah k. A similar estimate holds for the derivatives of P^ since -^P^{x) = 2XP^^^{x), 
see [24]. 

The set of orthonormal ultraspherical functions p'^{9) ~ sin^ 9P^{cos9) is a basis of L^{{0, tt), dx). 
The ultraspherical functions are eigenfunctions of the differential operator 

d^ A(A-l) 

namely, l\p^{9) = {k + \)'^p'^{9). By using the estimates for P^ given above, we can easily check 
that this system satisfies (5.1). Moreover, the heat-diffusion semigroup e^*'^ is positivity-preserving. 
This last assertion can be deduced directly from the facts that the heat-diffusion semigroup for the 
ultraspherical polynomials e~*^^ is positivity preserving, see [4], and e~*'^ = (C/o W) o (e~*^^) o (C7o 
W)~^, see Subsection 5.3.2 below. 

5.3.1. Proof of Harnack's inequality for {l\Y ■ We do this as we did for {T-LdY above by following 
parallel arguments. Let / G Dov[\{l\), / > 0, such that {l\Y f = in L,^{Ijd9), for some interval 
/ C (0,7r). Let u be the solution to the extension problem for l\ and this /. By the estimates 
mentioned above, u is classical. The potential here is V{9) — ^^^^g , which is a locally bounded 
function. Observe that -^p^{d) = —2\pj^'^\{9) + Xcot9p'^{9). Since cot is bounded in /, by following 

the same arguments as those for {'HdY i ^^ can get ||^'u(0,2/)||L2(7,de) — > ||./'(^)||L2(7.de), as y — > 0+. 
The conclusion follows by Theorem 3.2. 

5.3.2. Proof of Harnack's inequality for {L\Y ■ This is achieved by applying the transference method 
with M{9) = sin^ 9 and h{9) = 9. It readily follows that {LxY = {U o W)-^ o {l^Y o {U o W). 



6. Laplacian and Bessel operators 

In this section we will prove Theorem A for the fractional powers of the Bessel operator. This 
operator is a generalization of the radial Laplacian. For the sake of completeness and to show how the 
proof works, we present first the case of the fractional Laplacian on R", for which the more familiar 
Fourier transform applies. 

The main difference with respect to the examples given before is that these operators have a 
continuous spectrum and the Fourier and Hankel transforms come into play. 
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6.1. The Laplacian on M". Consider the fractional Laplacian defined by (— A)'^/(^) = |^| "^ /{£,), 

where f denotes the Fourier transform: f(f) = ceif) ~ - — ^ — -rr / f{x)e~"'^dx, f £ M". The 

* (27r)"/2 7r„ 

eigenfunctions of —A, indexed by the continuous parameter ^, are ^^{x) = e~*^'^, x e M", and 

{—A)ip^{x) — \£,\'^(fi^{x). Note that for any compact subset K C R" and any multi-index /? G Ng, 

\/3\ < 2, we have 



(6.1) 



\D^ 



n 



L-^iK) 



<ICI 



l/3| 



For any / e L^{K, dx), the heat semigroup is defined by e /(x) 



As 
(6.2) 



(27r)«/2 



-*l«l' 



cc(/)<P-c(a;)dC- 



e'^fix) <C 



-m- 



cdf)^.^^) rf^<Ci-"/'ll/IL.(K,d.), 2:eir, 



the integral that defines e^^f{x) is absolutely convergent in X x (0,T) with T > 0. Moreover, e*'^ is 
positivity-preserving in the sense of (2.2) because it is given by convolution with the Gauss- Weierstrass 
kernel. Note that, in this spectral language, Dom(-A) = {/ € i^(M",da;) : |CP/(0 e L'^{W,dx)} = 
{/ G L2(M",da;) : D^/ G i^^j^n^^^;)} = M^2,2^]gn)^ ^j^^ Sobolev space of functions in L2(M") with 
Hessian DV in ^^(R")^ 

Let us show Theorem A for (-A)'". Assume that / G W'^''^{W), / > and {-A^f = in 
L^{BR,dx), for some ball Br C M". By Theorem 3.2, we just must check that ||V2;ii(a-, y)j|i2(5^^2.) 
remains bounded as y — ?> 0"*". To that end, observe that for any x G Bj^ and y > 0, by (6.2), 

dt 










|e*''/(^K^lTTT^<C||/IL.(B„,,.)/ i 







^l+c 



< F{y), 



for some function F{y). This means that we can interchange integrals in u to get 



(6.3) 



u{x,y) 



y 



4'^r(cr)(27r) 



n/2 



ce (/)¥'-? (a^) 



_t|fH _«_ 



^i+a 



de 



By (6.1) and using the same arguments as above, it is easy to see that this double integral defines a 
function in C^{Bj^ x (0, oo)). So in this case w is a classical solution of (2.4). By using Plancherel's 
Theorem and (6.3) we have 

2 



\d,Mx,y) ~ fixMl. 



\dx) 



{d,^iuix,y)~fix))) iO 



LV 



(6.4) 



1 



(27r)'^ 

1 

(2^ 



H^Mif) 



4'^r(cr) 



e-*l«l e-^ 



dt 



2 



|(-*0)ce(/)^-e(a;)r 



y 



2ct 



4'^r(a) 7o 



(e-*l«l^ - l) 



rfC 



1 e -1* 



^i+a 



de 



Observe that the expression in square brackets above is uniformly bounded in y and it converges to 
when y -^ 0+. Moreover, as / G W^''^(R.''), \\{-i^j)cdf)\\mR'^,di) " W^xJWl^r^Mx) < oo. Hence, 
by dominated convergence in (6.4), dx u{x,y) converges to dx f in L'^{M.^,dx) as y ^ 0+. Whence 
Vxu{x, y) -J> '^xfi.x) as y -)► 0+, in L'^{Br, dx). 



6.2. The Bessel operators on (0,oo). Let A > 0. Let us denote by Aa the Bessel operator 

d^ 2A d 

Ax^--r^ 1~' a: > 0, 

dx"^ X dx 

which is positive and symmetric in L'^ {{0,oo), dm \{x)), where dm\{x) = x'^^dx, see [3, 18]. If 2A = 

n — 1, n G N; then we recover the radial Laplacian on R". Let Ji, denote the Bessel function of 

the first kind with order z/ and let us define (Pf{x) = x~^{^xy/'^Jx-i/2{^x), x,^ & (0,oo). Then, 

Axipf{x) = ^'^(p}{x), see [3]. These functions will play the role of the exponentials e~"^ in the case 

of the Laplacian. 
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We also consider the Bessel operator 

__il A' -A 

which is positive and symmetric in L^((0, oo), dx). Observe that the potential V{x) 

locally bounded function. If we let ipc{x) = x^(pc{x) then S\tl)f{x) = ^'^ip^{x), see [3]. The Hankel 

transform 

/•CJO 

/^ / i^\Cx)fix)dx 
Jo 
is a unitary transformation in L^((0, oo), da;), see [27, Chapter 8]. On the other hand, it is known 
that for any compact subset K C (0,oo) and k e Nq, there exist a nonncgative number s — SK,k 

and a constant C = CK,k such that \\ip^{x)\\L'^(^K,dx) < C, and |l^V'|(a;)llL~(K,da;) < C'|$|^ see [17]. 
Therefore parallel to the case of the Laplacian we can define the heat semigroup as 



e 



-ts^ 



f{x)= e-*«C5(/)^|(a;)de, 







where c^(/) = / f{x)ip^{x)dx. Moreover, 
Jo 

/■oo 

|e-*^V(-^)| < / e-'i\^{f)ij^{x) d^ < Ct-^l^ yWL-i^KA.) . ^ e ^' 

J 

so the integral that defines e~*^^f{x) is absolutely convergent in K x (0,r) with T > 0. Since e~*'^^ 
is positivity-preserving (sec [3] ) , we can follow step by step the arguments we gave for the case of the 
classical Laplacian to derive Theorem A for the operator {S\)°' . 

In order to get Theorem A for (A^)*^ we apply the transference method. Indeed, an obvious 
modification of Lemma 4.1 is applied with M{x) = x^ to get (Aa)*^ = U^"^ o {S\Y ° U. 

Acknowledgements. We are very grateful to Jose L. Torrea for many fruitful discussions. 

References 



[r 

[9: 
[lo: 

[11 

[12: 

[13: 
[14: 



I. Abu-Falahah, R. A. Maci'as, C. Segovia and J. L. Torrea, Transferring strong boundcdncss among Lagucrre 

orthogonal systems, Proc. Indian Acad. Sci. Math. Sci. 119 (2009), 203-220. 

I. Abu-Falahah and J. L. Torrea, Hermite function expansions versus Hermite polynomial expansions, Glasg. Math. 

J. 48 (2006), 203-215. 

J. Betancor, J. Dziubahski and J. L. Torrea, On Hardy spaces associated with Bessel operators, J. Anal. Math. 

107 (2009), 195-219. 

S. Bochner, Sturm-Liouville and heat equations whose eigenfunctions are ultraspherical polynomials or associated 

Bessel functions, in: Proceedings of the Conference on Differential Equations (dedicated to A. Weinstein), 23-48, 

University of Maryland Book Store, College Park, Md., 1956. 

L. Caffarelli, Some nonlinear problems involving non-local diffusions, in: ICIAM 07-6th International Congress 

on Industrial and Applied Mathematics, 43-56, Eur. Math. Soc, Ziirich, 2009. 

L. Caffarelli, S. Salsa and L. Silvestre, Regularity estimates for the solution and the free boundary of the obstacle 

problem for the fractional Laplacian, Invent. Math. 171 (2008), 425-461. 

L. Caffarelli and L. Silvestre, An extension problem related to the fractional Laplacian, Coram. Partial Differential 

Equations, 32 (2007), 1245-1260. 

L. Caffarelli and A. Vasseur, Drift diffusion equations with fractional diffusion and the quasi-gcostrophic equation, 

Ann. of Math. (2), 171 (2010), 1903-1930. 

E. B. Davies, Heat Kernels and Spectral Theory, Cambridge Tracts in Mathematics 92, Cambridge Univ. Press, 

Cambridge, 1989. 

E. B. Fabes, C. E. Kenig and R. P. Serapioni, The local regularity of solutions to degenerate elliptic equations, 

Comm. Partial Differential Equations, 7 (1982), 77—116. 

D. Gilbarg and N. S. Trudinger, Elliptic Partial Differential Equations of Second Order, Classics in Mathematics, 

Springer- Vcrlag, Berlin, 2001. 

C. E. Gutierrez, Harnack's inequality for degenerate Schrodinger operators, Trans. Amer. Math. Soc. 312 (1989), 

403-419. 

C. E. Gutierrez, On the Riesz transforms for Gaussian measures, J. Funct. Anal. 120 (1994), 107—134. 

C. E. Gutierrez, A. Incognito and J. L. Torrea, Riesz transforms, (^-functions, and multipliers for the Laguerre 

semigroup, Houston J. Math. 27 (2001), 579-592. 



HARNACK'S INEQUALITY 15 

[15] T. Lin, Y. Y. Li and J. Xiong, On a fractional Nironberg problem, part I: blow up analysis and compactness of 

solutions, preprint arXiv:1111.1332 (2011), 59 pp. 
[16] N. S. Landkof, Foundations of Modern Potential Theory, Die Grundlehren der mathematischen Wissenschaften, 

Band 180, Springer- Verlag, New York-Heidelberg, 1972. 
[17] N. N. Lebedev, Special Functions and Their Applications, Prentice-Hall, Inc., Englewood Cliffs, N.J., 1965. 
[18] B. Muckenhoupt and E. M. Stein, Classical expansions and their relation to conjugate harmonic functions. Trans. 

Amer. Math. Soc. 118 (1965), 17-92. 
[19] L. Roncal and P. R. Stinga, Fractional Laplacian on the torus, preprint (2012), 16 pp. 
[20] W. Rudin, Functional Analysis, McGraw-Hill Series in Higher Mathematics, McGraw-Hill Book Co., New York- 

Diisseldorf-Johannesburg, 1973. 
[21] M. F. Shlesinger, G. M. Zaslavsky and J. Klafter, Strange kinetics. Nature 363 (1993), 31-37. 
[22] L. Silvestre, Regularity of the obstacle problem for a fractional power of the Laplace operator, Comm. Pure Appl. 

Math. 60 (2007), 67-112. 
[23] P. R. Stinga and J. L. Torrea, Extension problem and Harnack's inequality for some fractional operators, Comm. 

Partial Differential Equations 35 (2010), 2092-2122. 
[24] G. Szegd, Orthogonal Polynomials, Fourth edition, American Mathematical Society Colloquium Publications 

XXIII, American Mathematical Society, Providence, R.I., 1975. 
[25] J. Tan and J. Xiong, A Harnack inequality for fractional Laplace equations with lower order terms. Discrete Contin. 

Dyn. Syst. 31 (2011), 975-983. 
[26] S. Thangavelu, Lectures on Hermite and Laguerre Expansions, Mathematical Notes 42, Princeton University Press, 

Princeton, NJ, 1993. 
[27] E. C. Titchmarsh, Intoduction to the Theory of Fourier Integrals, Third edition, Chelsea Publishing Co., New 

York, 1986. 
[28] N. S. Trudinger, Linear elliptic operators with measurable coefficients, Ann. Scuola Norm. Sup. Pisa (3) 27 (1973), 

265-308. 
[29] K. Yosida, Functional Analysis, Classics in Mathematics, Springer- Verlag, Berlin, 1995. 

Departamento de Matematicas y Computacion, Universidad de La Rioja, 26004 Logrono, Spain 
E-mail address: pablo-raul. stinga9uiiirioja.es 

School of Mathematics and Statistics, Wuhan University, 430072 Wuhan, China, -AND- Departamento 
DE Matematicas, Facultad de Ciencias, Universidad Autonoma de Madrid, 28049 Madrid, Spain 
E-mail address: zaoyangzhangchao@163.com 



